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Abstract 
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We study the influence of non-equilibrium phase space effects on the dynamics of 
CnJ ■ heavy ion reactions within the relativistic BUU approach. We use realistic Dirac- 

Brueckner-Hartree-Fock (DBHF) mean fields determined for two-Fermi-ellipsoid 
| configurations, i.e. for colliding nuclear matter, in a local phase space configura- 

tion approximation (LCA). We compare to DBHF mean fields in the local density 
approximation (LDA) and to the non-linear Walecka model. The results are further 
compared to flow data of the reaction Au on Au at 400 MeV per nucleon measured 
by the FOPI collaboration. We find that the DBHF fields reproduce the experiment 
if the configuration dependence is taken into account. This has also implications on 
i , the determination of the equation of state from heavy ion collisions. 
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1 Introduction 



A principal object of the investigation of heavy ion collisions at intermediate 
energies is to determine the nuclear equation of state (EOS), i.e. the proper- 
ties of nuclear matter in equilibrium away from saturation and at non-zero 
temperatures. However, the phase space distribution in a heavy ion collision 
is out of global and even local equilibrium through most of the collision time. 
Thus transport models have been developed to describe the evolution of the 
phase space [1—5]. Although these models are quite successful to reproduce 
data the attempts to determine the EOS, e.g., from the transverse flow in 
heavy ion collisions have not led to generally accepted results. Recently also 
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experimental collaborations, like FOPI or EOS, again focused on this problem 
[6,7]. 

The difficulty is in part due to the non-equilibrium effects mentioned above. 
The common practice of theoretical calculations is to apply phenomenological 
mean fields like Skyrme forces [2] in a non-relativistic or the Walecka model [8] 
and its non-linear extensions in a relativistic approach [4,9]. These forces con- 
tain parameters which allow to vary the characteristics of the corresponding 
EOS. In some cases a phenomenological momentum dependence taken from 
the empirical nucleon-nucleus optical potential has been added [2,10]. How- 
ever, non-equilibrium features of the phase space are not contained in these 
forces, a fact which certainly will lead to uncertainties when conclusions on 
the ground state EOS are drawn. 

A derivation of a kinetic equation from a microscopic non-equilibrium many- 
body theory shows that a consistent treatment requires to determine the ef- 
fective interaction in the nuclear medium for the non-equilibrium phase space 
configurations of heavy ion collisions [11]. Generally, the Dirac-Brueckner- 
Hartree-Fock (DBHF) theory provides a very successful approach to the many- 
body problem in nuclear matter [12,13] and in finite nuclei [14-16]. However, 
a solution of the Bethe-Salpeter equation for arbitrary anisotropic momentum 
configurations is extremely difficult. Approximate treatments should retain 
the most important features of such configurations. To improve on the local 
density approximation (LDA) the phase space has been represented locally 
by two Fermi ellipsoids separated by a relative momentum. We have called 
this the local configuration approximation (LCA) [17]. In contrast to the LDA 
which refers to ground state nuclear matter the LCA describes colliding nu- 
clear matter (CNM) and should be able to reproduce the main features of the 
anisotropic momentum space in a heavy ion collision over the entire reaction 
time. 

However, even for the idealized momentum configuration of two Fermi ellip- 
soids the Bethe-Salpeter equation has not yet been solved. Only in a non- 
relativistic framework this problem has been treated before [18]. Therefore 
as a first step Sehn et al. [19,20] developed a procedure to construct self- 
energies for colliding nuclear matter from an appropriate parametrization of 
DBHF ground state results. These self-energies take into account specific non- 
equilibrium features of the momentum space and approximately include the 
correlations of the T-matrix. Rather than introducing model parameters as, 
e.g. in the Skyrme or the Walecka model, the DBHF mean fields are connected 
in a parameterfree way to the free NN-interaction. In this sense we call these 
fields "realistic". 

We compare the present approach to experiment, i.e. to FOPI flow data of 
the system Au on Au at 400 A.MeV [6,21]. In addition we compare to the 
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LDA and to a standard force used in relativistic heavy ion calculations, the 
non-linear Walecka model (NL2). 



2 DBHF Mean Fields in the RBUU Approach 

As shown, e.g., in Ref. [11] the self-energy E = S s — 7^S^ which enters into 
the relativistic Boltzmann-Uehling-Uhlenbeck (RBUU) equation 

f(x,p*) = I c (1) 

via effective masses m* = M + ReS s , kinetic momenta = p^ + ReS^, and 
the field strength tensor F^ v in principle has to be determined from the corre- 
sponding T-matrix in non-equilibrium. Consistently with the kinetic equation 
(1) we represent £ in Hartree form by seperating off a linear density depen- 
dence through the definition of dynamical coupling functions 

iteE[/] = -T s [f] Ps [f] + ^r [f]f[f\ ■ (2) 

In Eq. (2) p s and are the scalar density and the baryonic four-current, 
repsectively, and [/] denotes the dependence on the phase space distribution 
/. The self-energy, Eq. (2), is of the same structure as the mean field in the 
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dcj-model [8], however, the coupling constants for the scalar (a) and 

the vector (uj) meson are replaced by the dynamical coupling functions r S) o[/]. 
These are Lorentz scalars and given by invariants of the T-matrix averaged 
over / [20]. In the LCA the r s , Eq. (2), are locally approximated by the 
values for the corresponding colliding nuclear matter configurations, and thus 
depend on the collective parameters of the CNM configuration, i.e. the Fermi 
momenta and the relative velocity of the local currents 

r S)0 [/] ' — ► r^ 2) (p Fl , p F2 , v Ic i) . (3) 

In this formulation the LDA means to use T(pp tot ) corresponding to a single 

— (12) 

Fermi sphere at the respective total density. The calculation of the T s is dis- 
cussed in detail in Ref. [20] . They are constructed by extrapolating the density 
and momentum dependent coupling functions T s (p,p B ) defined as in Eq. (2) 
for DBHF nuclear matter calculations to CNM configurations. There we used 
the DBHF results of Ref. [13] (without A-particles). Thus the T=0 EOS is 
the one given in Fig. 3.8 of Ref. [13], which gives a good fit of the saturation 
properties, an incompressibility of 250 MeV, and a good reproduction of the 
energy dependence of real part of the optical potential up to about 400 MeV. 



p^dl + (p*F"" + m*d%m*) d p * 
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In the LCA the collective parameters of the CNM configurations are deter- 
mined in every time step from the actual phase space distribution. To do so 
/ is decomposed into the respective contributions from projectile and target 
/ = + f( 2 \ The Fermi momenta p Fi = (2/37r 2 pg are defined in the rest 
frames of the currents by an invariant rest density [17,19] 



P^(x) = y/j^(x)j^(x) ,< = 1,2 (4) 



^%f { H^=uf{x)pf{x) . (5) 

The invariant relative velocity i> re i = \v re \\ is obtained from the streaming 
velocities u^} = (uq\u^) as 

. . . 4 2 'g"' - ^gw ... 



In the present calculations we assume symmetric configurations, i.e. p Fl = Vf 2 , 
which is a good approximation to the participant region. For highly asymmet- 
ric configurations where one density falls significantly below the other, what 
mainly occurs in the more peripheral reactions, the system is treated as one 
Fermi sphere in the LDA. 

To give an impression of the configuration effects we show in Fig. 1 the 
Schrodinger equivalent real part of the optical potential as a function of E\^. 
The DBHF nucleon optical potential in a nucleus-nucleus collision in the CNM 
approximation [20] is compared to the corresponding nucleon-nucleus optical 
potential [13] at the same density. In CNM the energy dependence of U opt 
originates from the dependence of the mean field on f rc i(-Eiab) whereas in the 
LDA, also shown in Fig. 1, the self-energy is taken at the corresponding den- 
sity and kept constant. In the latter case ReU opt scales linearly with energy as 
is the case, e.g., also in the Walecka model [4,8]. The configuration dependence 
(LCA) significantly softens the optical potential leading to a similar behavior 
as in nuclear matter when the full momentum dependence of the self-energy 
is included [13]. At two times saturation density the nucleus-nucleus optical 
potential is even softer. More exactly speaking we use "cold" CNM configura- 
tions to parametrize the phase space distribution, i.e. two sharp Fermi spheres. 
More consistently one could use two Fermi spheres of finite temperature. In a 
forthcoming work we analysed the phase space distribution in such a way and 
we find temperatures of less than about 30 MeV. In Ref. [13] it is found that 
ReUopt is only weakly changed at these temperatures and thus the assumption 
of zero temperature fields appears justified. 
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For the solution of the relativistic BUU equation we use the relativistic Landau- 
Vlasov method [5] which represents the phase space distribution by covariant 
gaussian testparticles in coordinate and momentum space. A fully consistent 
treatment of the kinetic equation in principle requires to treat the in-medium 
cross section on the same footing as the mean field [11]. Hence in Ref. [22] an 
in-medium cross section has been derived for CNM from the imaginary part 
of the DBHF self-energies. However, in this work we focus on the mean field 
and apply the standard Cugnon parametrization [23] for the cross section. As 
shown, e.g., in Ref. [25] the transverse flow does not react very sensitively to 
the cross section and we do not expect the present results to be significantly 
altered by the use of a consitently determined in-medium cross section. 



3 The System Au on Au at 400 A.MeV 



The present system is particularly well adapted to study non-equilibrium ef- 
fects. The bombarding energy is high enough to result in distinct two-Fermi- 
ellipsoid momentum configurations in the initial phase of the reaction and the 
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Fig. 1. Energy dependence of the DBHF optical potential. The solid lines repre- 
sent the nucleon-nucleus optical potential taken from ref. [13] at saturation density 
/9 sat (diamonds) and 2p sat (squares). The dashed lines represent the corresponding 
nucleon optical potential in a nucleus-nucleus collision determined in the colliding 
nuclear matter approximation at subsystem densities Pq 1 ^ + = p sa t (diamonds) 
and p^+Pq 2 -* = 2p sat (squares). The dotted lines refer to the DBHF nucleon-nucleus 
optical potential obtained in a simple LDA with no momentum dependence included 
at densities p sa t (lower curve) and 2p sat (upper curve). 



5 



densities reached in the compression phase are sufficiently high to test the 
density dependence of the EOS. Impact parameter selected high quality data 
of nuclear flow are provided by the FOPI collaboration [6]. Fig. 2 shows the 




Fig. 2. Transverse flow per nucleon in units of the nucleon mass for the reaction 
Au on Au at 400 A.MeV. Four ERAT bins running from semi-peripheral (top left) 
to central collisions (bottom right) collisions have been selected. Non-equilibrium 
DBHF mean fields were used (solid curves) as well as DBHF mean fields in the LDA 
(dashed curves). The dotted curves correspond to the non-linear Walecka model NL2 
and the diamonds represent the FOPI data of Ref . [6] . 

in-plane transverse flow per particle p x /A in units of the nucleon mass for 
four different ERAT bins. ERAT is defined as the ratio of transversal to lon- 
gitudinal kinetic energy deposited in the forward center-of-mass hemisphere 
[6] and provides a measure for the impact parameter bins which run from 
semi-peripheral {b ~ 6 fm, 0.30 < ERAT < 0.43) to central {ERAT > 0.72) 
collisions. For the analysis of the reaction we use the geometrical part of the 
FOPI filter [24] which takes into account acceptance cuts for different mass 
fragments. We therefore generated fragments after the collision (60 fm/c) by 
a phase space coalescence model. The influence of the filter on the observables 
is similar as in Ref. [6]. For the original experimental data the crossover of 
the flow from negative to positive values is slightly shifted to positive center- 
of-mass rapidity values. This behavior is quantitatively not well understood 
but is supposed to be due to recoil effects or double hits [24]. After applying 
the filter analysis the crossover of the theoretical curves is slightly shifted to 
negative rapidities. For a better comparison of the flow, i.e. the slope of p x /A, 
in Fig. 2 we readjusted the crossovers to p x /A = for theory and experiment. 

It is seen from Fig. 2 that the non-equilibrium DBHF mean fields (LCA) are 
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able to reasonably reproduce the data. The configuration dependence of the 
fields weakens their repulsion, as seen in Fig.l, resulting in less flow compared 
to a pure LDA treatment. This configuration effect on the flow is strongly 
increasing with centrality and reaches a magnitude comparable to different 
equations of state, i.e. about 30% of the total amount of flow. This can be 
understood by the increasing overlap of the nuclei since only the participant 
matter feels configuration effects. For the spectator matter our description 
always results in a LDA treatment. We also find that over the entire impact 
parameter range NL2 yields a too small transverse flow compared to the ex- 
periment. 

0.14 
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Fig. 3. Total transverse flow averaged over the foreward cms hemisphere in units of 
the nucleon mass for the same reaction as in Fig. 2. Non-equilibrium DBHF mean 
fields were used (squares) as well as DBHF mean fields in the LDA (triangles). The 
circles correspond to the non-linear Walecka model NL2 and the diamonds represent 
the FOPI data of Ref. [6]. 

In Fig. 3 we show the total amount of flow < p x (ERAT) > integrated over the 
forward center-of-mass hemisphere. Again, the LCA calculations are close to 
the data and also reproduce best the shape over the entire ERAT range. The 
larger deviation in the more peripheral collisions is probably due to the fact 
that the CNM mean fields are always determined for symmetric configurations 
(Pf 1 = Pf 2 )- For the very central collisions (high ERAT values) the LDA 
completely fails and overestimates the data by about a factor of 2 but becomes 
more reliable with increasing impact parameter. The NL2 model again strongly 
underestimates the flow. 
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Table 1 

Mean directed in-plane transverse flow normalized to the center-of-mass projectile 
momentum per nucleon. Calculations were performed with DBHF mean fields in 
non-equilibrium (LCA) as well as in the local density approximation (LDA) and 
with the non-linear Walecka model NL2. The experimental value is taken from Ref. 

M 



DBHF 


NL2 


EXP 


LCA LDA 






P x ir / p CM °- 187 °- 232 


0.103 


0.173 



The above results are summarized in Table 1. There the mean in-plane directed 
transverse flow normalized to the center-of-mass projectile momentum per 
nucleon is given. The experimental value [21] has been obtained for the PM4 
event class corresponding to semi-central (3 < b < 5) collisions. This quantity 
provides a good measure of the global repulsion of the model. The DBHF 
result is close to the experimental value when treated in the local configuration 
approximation (LCA). The local density approximation (LDA) overestimates 
the experiment by about 30% and the NL2 parameter set of the non-linear 
Walecka model underestimates the data by about 40%. 

To understand the influence of the mean field on the reaction dynamics we 
emphasize that in a relativistic approach the major part of the repulsion is due 
to Lorentz forces generated by the vector field (see Fig. 1). The repulsion and, 
correspondingly, the stiffness of the EOS in first order is given by the effective 
mass; the smaller the effective mass the more repulsive is the model [3,10]. 
Thus NL2 with a large effective mass (m*/m =0.83) is relatively soft whereas 
the DBHF EOS (m*/m =0.586) is rather stiff. The momentum dependence 
of the mean field which is taken into account in the LCA in an averaged form 
actually weakens this repulsion as also found in Ref. [10]. The small vector 
fields of NL2 lead to an underestimation of the flow. Furthermore, the density 
dependence of NL2 is too small to produce a sufficient amount of flow just 
by compression. Thus the use of parametrizations with a stiffer EOS but with 
the same m* will not significantly improve on this [4]. 

The same reaction as above has been analyzed in the QMD approach in Ref. [6] 
using forces based on the non-relativistic ground state BHF G-matrix [18,25]. 
There it was found that the strong and repulsive momentum dependence of 
the G-matrix is necessary to reproduce the data of the more peripheral col- 
lisions. However, in the more central collisions the non-relativistic G-matrix 
lacks a sufficient density dependence and results in a too small flow [6]. This 
drawback is resolved in a relativistic approach where the T-matrix provides 
the appropriate balance between density and momentum dependence. 
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4 Summary and Conclusions 



Realistic DBHF mean fields have been used in a detailed flow analysis of 
heavy ion collisions at intermediate energies. We have taken into account non- 
equilibrium features of the phase space in a local configuration approxima- 
tion using mean fields determined for a colliding nuclear matter scenario. The 
ground state DBHF self-energies are rather repulsive relative to, e.g., stan- 
dard parametrizations of the non-linear Walecka model. However, if the con- 
figuration effects are taken into account properly the DBHF mean fields are 
able to reproduce the overall behavior of the experimental flow data. Thus, 
the configuration effects turned out to be of major importance, especially if 
more central collisions are considered. A simple LDA treatment results in too 
repulsive mean fields. The intrinsic momentum dependence of the local con- 
figuration approximation weakens this repulsion considerably. An extension of 
the colliding nuclear matter approximation to anisotropic configurations and 
the inclusion of finite temperature effects and of an in-medium cross section 
may still improve the present results. A more detailed comparison to data 
including also longitudinal flow variables will be forthcoming. 

To summarize we conclude that DBHF mean fields based on realistic NN in- 
teractions are able to reproduce intermediate energy flow data, thus testing 
the DBHF approach also at high densities. It is found that the inclusion of 
non-equilibrium effects is important in order to learn something about the 
equation of state of ground state nuclear matter. 
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